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ABSTRACT 

The recently formulated Two Oscillator (TO) model interprets the lowest of 
the kilohertz frequencies of the twin peak quasi-periodic oscillations in X-ray 
binaries as the Keplerian frequency v^. The high twin frequency in this 
model holds the upper hybrid frequency relation to the rotational frequency of 
the neutron star's magnetosphere CI: v\ = v\ + 4(f2/27r) 2 . The vector f2 is 
assumed to have an angle 5 with the normal to the disk. The first oscillator 
in the TO model allows one to interpret the horizontal branch observed below 
100 Hz as the lower mode of the Keplerian oscillator under the influence of 
the Coriolis force, with frequency v L being dependent on u h , u K and 5. For 
some stars such as 4U 0614+09, Sco X-l and 4U 1702-42, v h , v K and u L have 
been observed simultaneously providing the opportunity to check the central 
prediction of the TO model: the constancy of 5 for a particular source. Given 
the considerable variation of each of these three frequencies, the existence of an 
observational invariant with a clear physical interpretation as a global parameter 
of the neutron star magnetosphere is an important test of the TO model. Using 
the results of recent observations of 4U 0614+09 we verify the existence of this 
invariant and determine the angle 8 = 15°. 6 ± 0.°5 for this star. The second 
oscillator in the model deals with a radial (presumably sound) oscillation and 
diffuse process in the viscous layer surrounding the neutron star. Our analysis of 
the viscous oscillation frequency vy and the break frequency of the diffusion 
shows that the spin value of the inner boundary of the transition layer for 4U 
0614+09 is at least two times more than values for 4U 1728-34 and Sco X-l. 

Subject headings: accretion, accretion disks — diffusion — stars:individual (4U 
0614+09, 4U 1728-34, Sco X-l, 4U 1702-42)— stars: neutron— X-ray:star— waves 
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1. 



Introduction 



This Letter contains further verification of the Two Oscillator (TO) model for kHz 
QPO's in NS binaries and the related low frequency features using the observational results 
for 4U 0614+09 from van Straaten et al. (2000) (hereafter VS00). In a series of the papers 
(Osherovich & Titarchuk 1999a, Titarchuk & Osherovich 1999, Osherovich & Titarchuk 
1999b, Titarchuk, Osherovich & Kuznetsov 1999, hereafter OT99a, T099, OT99b and 
TOK respectively) the authors offered the TO model which explained the variation of the 
frequency difference between two kHz QPOs suggesting that the upper one, namely v h is an 
upper hybrid frequency of the Keplerian oscillator under the influence of the Coriolis force 
and the lower kHz QPO is the Keplerian frequency uk- This identification is a principal 
difference between the TO model and other QPO models presented so far (e.g. Stella & 
Vietri 1998 Miller, Lamb, Cook 1998, Kaaret, Ford & Chen 1997) since in other models, 
the highest kHz QPO is identified as a Keplerian frequency. 

According to the TO model, the Keplerian oscillator has two branches characterized by 
high frequency Uh (~ 1 kHz) and by low frequency vl (~ 50 — 100 Hz). The frequency vl 
depends strongly on the angle S between the normal to the neutron star disk and Q— the 
angular velocity of the magnetosphere surrounding the neutron star. In the lower part 
of the QPO spectrum (~ 10 Hz), the second oscillator of the TO model describes the 
physics of the viscous transition layer, namely, radial viscous oscillations with frequency 
u v (previously called "extra noise component") and the diffusive process in the transition 
region (the innermost part of the disk) which is characterized by the break frequency 
Ub. According to the TO model, all frequencies (namely Vh, vl, and uy) have specific 
dependences on u K . Correlations of v h and u v with kHz frequencies are fitted in T099 by 
the theoretical curves using the dimensionless parameter, ax' 



a K = m(:ro/3) 3/2 (z/o/363 Hz) 



(1) 
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where normalized mass m = M/M e , radius of the inner edge of disk x = Rq/Rs m 
Schwarzchild radius (i?s units) and u Q is the spin of the disk inner edge. 

Van Straaten et al. (2000) have found for 4U 0614+09 that dependences of u v and u b 
on v K do not follow the theoretical curves which fit the data points for 4U1728-34 (T099) 
and Sco X-l (TOK). In this Letter, we revisit this issue and show that the observational 
correlations for source 4U 0614+09 are related to the theoretical ones for the same 
parameter a K = 1.03 (as for 4U 1728-34 and Sco X-l) but with a spin of the inner edge of 
accretion disk higher than 363 Hz. Such a high spin leads to the high compactness of the 
central object because 



We check predictions of the TO model regarding vy vs z/k dependence as well as relation 
between u b and u K for all available observations (see §3). 

In OT99b, the authors demonstrated for the source 4U 1702-42 that the inferred angle 
5, (see Eq 2, 5 there) 



depends only, on z/k, vk an d vl and stays the same (3.9° ± 0.2°) over significant range of z/k 
(650-900 Hz). In Formula (3), 5 is assumed to be small. A general formulation for any 5 is 
presented, in OT99a with conjecture that eventually 5 as an invariant will be found for all 
20 sources known to have QPOs. Using series of observations where Vh, v K , ul are detected 
simultaneously (van der Klis et al. 1997; Markwardt, Strohmayer & Swank 1999; van 
Straaten et al. 2000) we check that the angle 5 is a true invariant for three specific sources: 
Sco X-l, 4U 1702-42, 4U 0614+09 (see §2). We will also present the magnetospheric 
rotational profile for 4U 0614+09 inferred from VS00 data in §4. Discussion and summary 
follow in the last section. 





5 = arcsin {v\-v^) 1/2 (u L ^h/^K) 



(3) 



2. 5— Invariant and Verification of Two-Oscillator Model 
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We have used the frequencies measured for sources 4U 0614+09, Sco X-l and 4U 
1702-42 where the low and high kHz QPO peaks u K and v h and HBO frequencies v L are 
measured simultaneously. The resulting values of 5 calculated from Eq. (3) are shown in 
Figure 1. Indeed, for each of these sources, the 5— values show little variation with z/ K , 
v h , v L . They are 3°.9 ± 0.2 and 5°.5 ± 0.5 for 4U 1702-42 and Sco X-l respectively. For 
the source 4U 0614+09 the point at z/k — 821.6 Hz is not included since the z/k peak is 
very broad with FWHM= 205 Hz, which is almost order of magnitude larger than that for 
any other points for the same star according to Table 1 in VS00. The angle 5 obtained for 
the source 4U 0614+09 is 15°. 6 ± 0°.5 and is significantly larger than for other stars. The 
existence of the 5— invariant predicted by the TO model (OT99b) is a challenge for any 
other QPO model. It is important to note that all frequencies included in this 5— relation 
are observed frequencies and thus the relation Eq. (3) is a model independent invariant. 

3. Break Frequency vs kHz QPO Frequency Correlations 

Further tests of the TO model can be done using a comparison of the observed 
correlation of break and kHz QPO frequencies (VS00) with the theoretical dependences 
derived in Titarchuk, Lapidus & Muslimov (1998, hereafter TLM) and T099. In Fig. (2) 
(upper panel), we present the theoretical curves calculated using Eq. (9) in T099 for 
different values of ^0,363 — ^o/363 Hz. For the source 4U 0614+09 we found the Lebesgue's 
measure (for definition, see TOK) is better than 15% for ax = 1.03 and ^0,363 > 2 which is 
at least 3 times less than that for z/ ,363 = 1 (see upper panel in Fig. 2). We fitted the data 
points by the curves with z/o,363 > 3 [uq > 1089 Hz) and we have found that the quality of 
those fits is improved only slightly. Thus it is difficult to get the constraints on vq better 
than u > 800 Hz. Recently Bhattacharyya et al. (2000) using the full General Relativity 
technique argued that the NS in Cyg X-2 rapidly rotates with the frequencies of order of 1 
kHz which is very close to that we obtained here for 4U 0614+09. 
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We approximate the solution of Eq (9) of T099 by a polynomial for a K = 1-03 and 
^0,363 = 1, namely, 



where B 1 = -1.68 x KT 4 , B 2 = 1.36 x KT 5 Hz" 1 , B 3 = -2.47 x 1(T 8 Hz" 2 , B A = 4.45 x KT 11 
Hz -3 . We fit the observational break frequencies by the theoretical curve which is a function 
of parameter z/o,363- For a given z/o,363 that curve can be obtained by the scaling the 
argument x of the polynomial P±{x): 



The normalization of the curve is controlled by the constant C b which reflects the properties 
of the specific source. For 4U 1728-34 the observed frequencies are fitted by the curve with 
a K = 1.03 (T099) for which the constant C b — \ with R — 11 km for 1.4 solar masses. For 
4U 0614+09 we fit the data points by the curves with the parameters ^0,363 = 1, 2, 3, (i.e. 
v>o = 363 Hz, 726 Hz and 1089 Hz) for which the corresponding constants C& = 2, 22.8, 55.4 
are chosen respectively. Using Eq (2), one can estimate mass-radius (MR) relations for a 
given value of ^0,363- For u 0j363 = 2 and 3 the radius R ~ 6.4 and 4.9 km respectively, 
for the standard NS mass, 1.4 solar masses. There is definitely an indication of higher 
compactness for 4U 0614+09. It is true that there is a possibility to elaborate a much more 
stringent condition on the MR relation using Eq. (2) and fits of the data in Fig. 2 which 
may lead to small masses and radii for the star — m < 1 and R < 9 km. In that case 
such a low mass compact star can be easily explained in terms of strange star models (e.g. 
Bombaci 1997). In fact, under constraints imposed by the TO model Li et al. 1999 pointed 
out the possibility that 4U 1728-34 is a strange star rather than a neutron star. Thus our 
arguments of the high compactness in 4U 0614+09 support a strange star concept. But 
the final conclusion on MR constraints cannot be worked out without the full relativistic 
treatment which is out of the scope of this paper. 



v* = C b ■ P A (u K ) = C 6 (5iz/ K + B 2 v 2 K + B 3 ul + £ 4 *4) 



(4) 



Vbfa) = C bP*(vK/ '^0,363)- 



(5) 
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According to the low panel in Fig. 2, v h is not lower than 400 Hz for any value of is K . 
Because 

^=[4 + 4(^/27r) 2 ] 1/2 , (6) 

the rotational magnetospheric frequencies v mag = (fl/2n) should have an asymptote for 
large radii (or small values z/r ^ 2v mag ) which equals approximately 200 Hz. The detailed 
profile of u mag constructed below confirms this expectation. 

4. Inferred Rotational Frequency Profile of the NS Magnetosphere 

From the observed kHz frequencies (VS00, Table 1), v h and z/ K , the profile of 
v mag = Q(v K )/2n has been calculated according to formula (6) and modeled using 
theoretically inferred magnetic multipole structure of differentially rotating magnetosphere 
(OT99a) 

Vmag = C + C X vf + + C 3 4 (7) 

where C 2 = 2(C 1 C 3 ) 1 / 2 . The constants C = v G mag = 190 Hz, C x = 6.52 x 10~ 2 Hz" 1 ^ 
C 2 = -1.19 x 10~ 5 Hz" 5 / 3 and C 3 = 0.91 x 10" 9 Hz" 3 have been obtained by the 
least-squares fit with x 2 — 29.8/25 (OT99a). The resulting fit is shown in the lower panel 
of Fig. 3 by solid line. The v mag profiles for Sco X-l (which is very similar to those for 
4U 1608-52 and 4U 1728-34) and that for 4U 0614+09 differ in the sign of curvature: it is 
negative for the profile of Sco X-l profile and positive for that of 4U 0614+09. The v mag 
profiles for Sco X-l, 4U 1608-52 (as well as for 4U 1728-34) according to analytical formula 
(14) in OT99a have a maximum at large z/k (small radii) and also have a minimum at small 
z/k (large radii). The adjustment of the Keplerian disk to the sub-Keplerian rotation of the 
NS in 4U 0614+09 occurs at much larger distances from the central source (3.6-7 NS radii) 
than that for 4U 1728-34 where the adjustment radii are within 1.6-3.4 NS radii. Thus the 
blobs (formed at the adjustment radius as result of the super— Keplerian motion [TLM]) 
presumably probe the different parts of the magnetosphere: in the 4U 0614+09 case, it is 
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the outer part of the magnetosphere, near the minimum of the rotational frequency while 
in the 4U 1728-34 case, it is a part which is closer to the inner part, near the maximum. 
Thus one should make an expansion of u mag over the Chandrasekhar potential A near the 
minimum in the former case and near the maximum in the latter case (see OT99a for 
details). We may try to predict v mag for much smaller u K (around 150 Hz) than those 
frequencies z/k (between 400 and 850 Hz) where we have constructed the profile (7). Using 
VS00 data points for 8 August with u K = 153.4 ± 5.6 Hz and v h = 449.4 ± 19.5 Hz and Eq. 
(6), we find u rnag = 211. ± 23. Hz which is in good agreement with the value v mag = 235 Hz 
obtained from the analytical expression, Eq. (7). 

As shown in Fig. 3 (the upper panel) the difference in v mag profiles for Sco X-l and 
for 4U 0614+09 results in different profiles for Av = Uh — i^k- namely, Av monotonically 
decreases with v K for Sco X-l but has a minimum for 4U 0614+09. According to VS00 
Av = 311.8 ± 1.8 Hz except one point (not included in Fig. 3) which deviates 11 a from 
this trend. We will return to this point in another paper. 

5. Break Frequency vs Low Lorentzian QPO Frequency Correlation 

According to solutions presented in T099 and TOK, the break frequency v b correlates 
with the low Lorentzian frequency vy (which we interpreted as a viscous frequency in the 
TO model) through power law relation 

v b « 0.041z/£ 61 . (8) 

It was confirmed in TOK using the data of Wijnands & van der Klis (1999) that the same 
index approximately is valid for the wide class of NS and BH objects. In Figure 4 we 
demonstrate that the observed correlation of v b vs v v in 4U 0614+09 is described well by 
the same power law dependence and the quality of this fit is similar to that of 4U 1728-34, 
i.e. the Lebesgue's measure is less than 20%. Equations (5) and (8) give us theoretical 
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dependences of v v on is K , assuming z/o,363 = 3. For August 8, 1996 the theoretical value 
v v w 9 Hz with inherited 20% accuracy comes close to the observed vy = 11.2 ± 1.4 Hz as 
given in VS00. 

6. Discussion and Conclusions 

In this Letter we have verified predictions of the two-oscillator (TO) model for three 
sources. The unifying characteristic of spectra for both oscillators which share the common 
boundary at the outer edge of the viscous transition layer is the strong dependence of all 
frequencies in the model on v^. The angle 5 as a global parameter describes the inclination 
of the magnetospheric equator to the equatorial plane of the disk. Measured locally for 
different radial distances (therefore different i/k), S may vary considerably if the observed 
oscillations do not correspond to the predicted low Keplerian branch. The constancy of 6 
shown in Figure 1 argues in favor of the TO model. The parameter 5 is critical in evaluating 
the differences between the spectra of different sources and we suggest that the higher v L 
for 4U 0614+09 and 4U 1728-34 (than that for Sco X-l) are mainly a result of a larger 
angles 5. 

We consider the angle 8 to be a fundamental geometric parameter for neutron stars 
and may be related to the evolution and prehistory of these systems. During the NS life 
time, the NS is believed to spin up gaining more angular momentum from the companion 
and the axis of magnetosphere may change with respect the rotational axis of the star (or 
the disk) (Bhattacharya 1995). The X-ray bursters are found in old population (age of 
order 10 9 years) concentrated towards the galactic center (Lewin, van Paradijs & Taam 
1995) and thus it is not surprising that in these systems the angle S between spin axis and 
the magnetic axis is small presumably as a result of long evolution. The absence of X-ray 
pulsations from the LMXBs and the fact that type I X-ray bursters never occur in systems 
which show pulsations (Lewin, van Paradijs & Taam 1995) probably can be explained by 
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the small angles 5. In contrast, in X-ray pulsars which are mostly observed in the High 
Mass Binaries (HMXBs), 5 is not small: the mean value of 5 is 27° (Leahy 1991) and the 
age of these systems is much smaller, of order of 10 7 years. 

The authors acknowledge discussions with Joe Fainberg, Sergey Kuznetsov and Craig 
Markwardt. We are grateful to Steve van Straaten for the data which enabled us to make 
comparisons with the data. We are also grateful the referee for the fruitful suggestions. 
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Fig. 1.— 5— invariant. Inferred angles 5 between the rotational axis of magnetosphere 
and the normal to the plane of Keplerian oscillations as a function of low peak kHz QPO 
frequency for three sources: 4U 0614+09, Sco X-l, and 4U 1702-42. The 5— invariant (eq. [3]) 
is calculated using the observational frequencies Vh, and vl only. 

Fig. 2. — Low panel: Break frequency of a broken power law function of the power spectrum 
density versus the frequency of the kHz QPOs for 4U 0614+09 (Straaten et al. 2000). 
Open circles are for the dependence versus low kHz peak and open squares are the 
v b — dependences versus high kHz peak. The solid line is the theoretical curve for break 
frequencies calculated using (T099 eqs. [4-5, 9, 13] corresponding to the parameter at = 1.03 
and vq = 1089 Hz. Upper panel: The theoretical dependences of the break frequency versus 
the low peak kHz QPOs for different values of the spin frequency of the inner edge of accretion 
disk vq. Observational points for 4U 0614+09 (Straaten et al. 2000) open circles and for 4U 
1728-34 (Ford & van der Klis 1998) solid circles 

Fig. 3. — Low panel: Inferred rotational frequency of the NS magnetosphere (solid line) as a 
function of the frequency of the low kHz QPO peak (Eq. 7 for 4U 0614+09 and see OT99a 
for Sco X-l). For 4U 0614+09 (van Straaten et al. 2000) observational points shown as 
solid circles and for Sco X-l (van der Klis et al. 1997) as open circles. Upper panel: The 
peak difference of kHz QPO frequencies for the same sources. The solid lines are theoretical 
curves obtained using the above theoretical dependences of the rotational frequencies. 

Fig. 4. — Break frequency versus main viscous frequency for 4U 0614+09. Solid line is the 
theoretical curve which is almost a power law with index 1.61. 
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